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Abstract. We present preliminary results from surface 
photometry of a sample of early-type galaxies observed 
with NICMOS on board HST in the F160W (i/-band) 
filter. Dust is found to occur mostly in galaxies which 
possess unresolved nuclei, and in a few cases the pres- 
ence of dust causes the isophotes to appear boxy. We per- 
formed a 2-D modeling of the host galaxy light distribu- 
tion using the empirical Nuker function, which proves to 
be an adequate representation of the surface brightness 
profile for most early-type galaxies. The accurate mod- 
eling of the bulge light allows us to obtain the magni- 
tude of the unresolved central source without contamina- 
tion from the underlying stellar population. Our classifi- 
cation of early-type galaxies into core-type and power- 
law galaxies agrees with the results from previous work 
based on HST optical images, with core profiles occurring 
in luminous galaxies and power-law profiles occurring in 
low-luminosity galaxies with disky isophotes. We do not 
find a dichotomy in the distribution of inner slopes. Un- 
resolved nuclear point sources are evident in about one- 
third of the sample galaxies, with magnitudes in the range 
13.5 <mH < 17.3. Most of the type 1 Seyfert and LINER 
nuclei (83%) reveal a distinct point source, while the de- 
tection rate is much less for type 2 nuclei (35%). The oc- 
currence of point sources appears to be independent of 
whether the underlying galaxy has core-type profile or 
power-law profile. 



1. Introduction 

Images of the central regions of galaxies at HST resolu- 
tion have made it possible to probe structures very close 
to the nucleus and to reveal unresolved point sources or 
compact star clusters. Optical images obtained using HST 
showed that the King model and the de Vaucouleurs r^/'' 
law which were conventionally used to describe the sur- 
face brightness profiles of early-type galaxies do not pro- 
vide a good fit for the regions within ~ 10" (Ferrarese 
et al. 1994; Lauer et al. 1995; CaroUo et al. 1997; Faber 
et al. 1997). The surface brightness profiles in the central 



regions can be parameterized by the empirical "Nuker" 
function, which is essentially a double power-law of the 
form 
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where (3 is the slope outside the break radius r;,, 7 is the 
inner slope, a controls the sharpness of the transition be- 
tween the outer and inner slopes, and h is the surface 
brightness at the break radius (Lauer et al. 1995). Early- 
type galaxies can be classified into core-type or power-law 
systems based on their surface brightness profiles. Core- 
type galaxies show a significant change in the slope of the 
surface brightness profile within the break radius rt with 
the inner slope 7 being much shallower than for the outer 
regions. Power-law galaxies have profiles that show no sig- 
nificant change in slope within the inner 10" and remain 
as steep power laws all the way to the resolution limit 
(r = O'.'l). A value of 7 < 0.3 is adopted as the criterion 
for classifying a galaxy as a core-type. Faber et al. (1997) 
found that no galaxies had 7 in the range 0.3 to 0.5 and 
therefore that there was a clear dichotomy in the distri- 
bution of 7 values between core galaxies and power-law 
galaxies. Interestingly, the surface brightness profiles de- 
rived from WFPC2 y-band images of early-type galaxies 
with kinematically distinct cores yielded values of 7 be- 
tween 0.3 and 0.5, thereby removing the dichotomy (Car- 
oUo et al. 1997). Surface photometry of early-type galaxies 
based on ff-band images from HST have been presented 
by Quillen, Bower, & Stritzinger (2000), along with the 
analytic Nuker fits, but they did not decouple the contri- 
butions from the bulge and point sources. 

In this work, we aim to model the galaxy surface 
brightness closely and to extract the point source mag- 
nitudes free from contamination by the bulge light, for a 
sample of galaxies that host AGNs. Parameterizing the 
bulge and nucleus as separate components enables us to 
probe the relationship between the properties of the host 
galaxy and the nuclear activity. The NIR images are most 
suitable for deriving the surface brightness profiles because 
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Table 1. Properties of sample galaxies^ 



Galaxy Hubble Type Bt Spectral Class 




NGC 221 


E2 


9.0 


-15.5 


A 


NGC 404 


SO 


11.2 


-15.9 


L2 


NGC 474 


SO 


12.3 


-20.5 


L2 :: 


NGC 524 


SO 


11.3 


-21.3 


T2 :: 


NGC 821 


E6? 


11.6 


-20.1 


A 


NGC 1052 


E4 


11.4 


-19.9 


L1.9 


NGC 2685 


SBO 


12.1 


-19.2 


S2/T2 


NGC 3115 


SO 


9.8 


-19.3 


A 


NGC 3379 


El 


10.2 


-19.3 


L2/T2 :: 


NGC 3384 


SBO 


10.8 


-18.7 


A 


NGC 3593 


so/a 


11.8 


-17.2 


H 


NGC 3900 


SO 


12.2 


-20.1 


L2 : 


NGC 4026 


SO 


11.6 


-19.5 


A 


NGC 4111 


SO spin 


11.6 


-19.5 


L2 


NGC 4143 


SABO 


11.6 


-19.2 


L1.9 


NGC 4150 


SO 


12.4 


-17.5 


T2 


NGC 4261 


E2 


11.4 


-21.3 


L2 


NGC 4278 


El 


11.0 


-18.9 


L1.9 


NGC 4291 


E 


12.4 


-20.0 


A 


NGC 4374 


El 


10.0 


-21.1 


L2 


NGC 4406 


E3 


9.8 


-21.3 


A 


i\Tr^ /I /1 1 7 
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brSU : spin 
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A 


NGC 4472 


E2 


9.3 


-21.8 


S2 :: 


NGC 4589 


E2 


11.6 


-20.7 


L2 


NGC 4636 


EO 


10.4 


-20.7 


L1.9 


NGC 5273 


SO 


12.4 


-19.2 


S1.5 


NGC 5548 


so/a 


13.3 


-21.3 


S1.5 


NGC 5838 


SO 


11.9 


-20.5 


T2 :: 


NGC 5982 


E3 


12.0 


-20.8 


L2 :: 


NGC 6340 


so/a 


11.8 


-20.0 


L2 


NGC 7457 


SO 


12.0 


-18.6 


A 


NGC 7626 


E ; pec 


12.1 


-21.2 


L2 :: 


NGC 7743 


SBO 


12.3 


-19.7 


S2 



^ Galaxy properties taken from Ho et al. (1997) 



of the reduced contamination from dust and young star 

clusters, and tlicy sample the underlying smooth stellar 
population better than optical images. 

2. The Sample and Data Analysis 

Early-type (E and SO) galaxies were selected from the sam- 
ple of galaxies (Table 1) studied spectroscopically by Ho, 
Filippenko, & Sargent (1997). Archival HSTmCUOS im- 
ages are available for 33 of these galaxies (14 E, 16 SO, and 
3 SO/a) from different snapshot survey programs. The data 
were obtained in the F160W filter using NIC2 and NIC3 
cameras. The NIC2 images have a field of view of 19'.'2 x 
19'.'2 with a scale of O'.'076/pixel. The NIC3 images have 
a field of view of 51'.'2 x 51'.'2 with a scale of O'.'2/pixel. 
The CalnicA images were used after performing additional 
steps suggested in the NICMOS handbook which include 
pedestal removal, creating bad pixel masks, and interpo- 
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Fig. 1. The results of isophotal analysis are shown for 
NGC 4374 (top) and NGC 3115 {bottom). The left image 
is the observed _ff-band image and the right image is the 
residual after subtracting the model derived from ellipti- 
cal isophotes. Each image is 10" x 10" and is centered on 
the galaxy. The orientation is arbitrary. The lower panels 
show the variation of ellipticity, position angle, and the 
shape parameter (B4) along the semi-major axis. 
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Fig. 2. Nuker-law fits to the core-type galaxy NGC 3379 
(top) and the power-law galaxy NGC 7626 {bottom). 
The solid points represent the observed surface brightness 
profile. The solid line is the 1-D profile extracted from 
isophote fits to the generated model galaxy. Upper limits 
to nuclear magnitudes are given. 

lating over the bad pixels (Dickinson 1999). The ramp 

fitting procedure used in the CalnicA removes most of the 
cosmic rays from the images. Residual cosmic rays, defec- 
tive pixels, and the coronographic hole (in NIC2 images) 
were masked during further analysis. 

3. Isophotal Analysis 

The ellipse task available in STSDAS was used to fit el- 
liptical isophotes to the galaxy images and thereby obtain 
the surface brightness profiles and the radial variation of 
ellipticity and position angle. The results of the isophotal 
analysis were used to create smooth model images of the 
galaxies, and residual images which highlight the high- 
spatial frequency features were obtained by subtracting 



Fig. 3. Nuker-law fits with point source for the core-type 
galaxy NGC 4374 (top) and the power-law galaxy NGC 
474 {bottom). The solid points represent the observed sur- 
face brightness profile. The solid line is the 1-D profile ex- 
tracted from isophote fits to the generated model galaxy. 
The dotted line is the intrinsic bulge profile obtained by 
omitting the point source in the model. 



the model images from the observed images. The presence 
of features such as dust and weak disks causes variations 
in ellipticity and P.A. The deviation from perfect ellip- 
tical isophotes can be expressed as a Fourier series, and 
the Fourier coefficient of the fourth cosine term in the 
expansion scries (B4) is negative for boxy isophotes and 
positive for disky isophotes (Jedrzejewski 1987). The re- 
sults from the isophotal analysis are shown in Figure 1 
for NGC 4374, in which dust is revealed in the residual 
image, and for NGC 3115, which has a dominant stellar 
disk. Notice that the isophotes tend to get boxy around 
r = 0'.'2 - C'S in the case of NGC 4374 due to the dust 
lanes. 
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4. 2— D Modeling Using Analytic Functions 

Earlier HST studies of the central surface brightness pro- 
files of galaxies involved deriving the central parameters 
by fitting the 1-D profiles with analytic functions. The un- 
ccrtainities caused by the PSF were overcome by perform- 
ing the fits on profiles derived from deconvolved images 
(Lauer et al. 1995) or by convolving the analytic 1-D fits 
with the PSF and then comparing to the observed profiles 
(CaroUo & StiaveUi 1998). Most of the data used in the 
present work were obtained through snapshot survey pro- 
grams and do not have the high S/N required to perform 
reliable deconvolution of the images. Since a more reliable 
approach would be to sample the PSF in two dimensions, 
we choose to parameterize the galaxy bulge through a 2- 
D modeling using analytic functions and convolving the 
model image with the 2-D PSF image. This was done using 
the GALFIT program (Peng et al., in preparation), which 
enables us to create model images for different analytic 
functions (e.g., Sersic/de Vaucouleurs, Nuker, exponen- 
tial, Gaussian) and convolve them with the PSF. In lieu 
of observing PSFs simultaneously, TINYTIM simulations 
provide a good alternative (Krist & Hook 1997). GALFIT 
allows modeling the contribution from unresolved point 
sources or resolved star clusters in addition to the un- 
derlying galaxy. The program outputs the fit parameters 
for the galaxy and the point source or star cluster, along 
with images of the model galaxy and residuals. Following 
the results from earlier HST studies of early-type galax- 
ies, we modeled the galaxy light using the empirical Nuker 
function. When appropriate, we also fit a compact source 
simultaneously in order to deblend a component in excess 
of the underlying galaxy component. When no compact 
source was required, we obtained an upper limit to the 
point source magnitude. A comparison of the 1-D surface 
brightness profiles derived by doing ellipse-fits on the ob- 
served image and the model image created by GALFIT, 
shows that the 2-D modeling reproduces the observed pro- 
files very well (Figures 2 and 3). 
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Fig. 4. The FP relation between ^f, and for core galax- 
ies. 

for their sample of ellipticals with kinematically distinct 
cores. Core galaxies obey the fundamental-plane (FP) re- 
lations in the (log r^, /ib)-plane with brighter galaxies hav- 
ing larger cores with lower surface brightness (Figure 4). 

5.2. Properties of Unresolved Nuclei 

About one-third of our sample galaxies have clear evidence 
for the presence of an unresolved central compact source, 
and the frequency of occurrence of point sources is al- 
most equal for core galaxies and power-law galaxies. How- 
ever, all the central sources associated with core galaxies 
are AGNs by their nuclear spectral classification, while 
power-law galaxies mostly contain absorption-line nuclei, 
transition objects, and weak AGNs. Most type 1 Seyfert 
and LINER nuclei (5 out of 6) have point sources com- 
pared to type 2 nuclei (5 out of 14). The unresolved nuclei 
have magnitudes in the range 13.5 < niR < 17.3, with the 
exceptionally bright Seyfert nucleus in NGC 5548 having 
niH = 11.9 mag. 



5. Results 

5.1. Properties of the Host Galaxies 

Early-type galaxies appear to have a smooth light distri- 
bution in the H band, showing dust features only in very 
few cases, mostly associated with galaxies hosting unre- 
solved nuclei (e.g., NGC 1052, 4150, and 4374). The pres- 
ence of dust in some cases is found to cause boxiness of 
the isophotes resulting in negative B4 values. We confirm 
the results from earlier studies that core galaxies are lu- 
minous elliptical galaxies with boxy or neutral isophotes, 
while power-law galaxies have low luminosities and disky 
isophotes (Faber et al. 1997). However, we do not find 
a dichotomy in the inner slope values (Ravindranath et 
al. 2001). Instead, we find a continuous distribution of 
7 values similar to that found by CaroUo et al (1997) 
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